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Problems  of  constructing  a  sediment  budget  for  an  inlet  is  often  a  challenge  due  to  constraints  on  data  collection, 
quality  of  data,  and  assumptions  one  is  willing  to  consider.  Existing  literature  does  not  provide  a  systematic 
foundation  for  equations  considered  and  often  engineers  and  geologists  typically  do  not  properly  consider  the 
mathematical  constraints  placed  on  the  problem  of  sediment  budget  construction.  The  present  technical  note 
presents  a  simple  methodology  for  assessing  sediment  budget,  i.e.  what  is  possible,  and  what  is  not  possible 
via  a  matrix  equation  system  and  linear  algebra.  It  is  hoped  that  engineers  and  geologists  will  not  impose 
unrealistic  expectations  on  the  sediment  budget  system  via  using  such  an  approach.  Some  simple  examples 
are  given  of  both  pitfalls  and  correct  approaches  to  sediment  budget  development. 

©  2011  Elsevier  B.V.  All  rights  reserved. 


1.  Sediment  budget  discussion 

The  problems  of  constructing  sediment  budgets  for  inlets  is  not 
new,  but  is  always  a  challenge  due  to  constraints  on  data  collection, 
quality  of  data,  and  assumptions  the  engineer  or  geologist  is  willing 
to  make.  Various  literature  exists  on  sediment  budget  construction 
(see  for  example,  Jarrett  (1991),  Bodge  (1999),  Rosati  and  Kraus 
(1999),  Dean  and  Dalrymple  (2002),  U.S.A.C.E.  (2003),  Rosati 
(2005)  and  others)  although  previous  discussions  have  not  typically 
addressed  the  mathematical  constraints  placed  on  the  problem  of 
properly  constructing  a  sediment  budget.  The  present  technical  note 
presents  another  way  of  looking  at  what  is  possible  and  what  is  not 
possible  in  sediment  budget  development  via  the  matrix  equations 
utilized  in  developing  the  sediment  budget  with  the  hope  that  engi¬ 
neers  and  geologists  will  be  careful  in  addressing  sediment  budget 
development  without  imposing  unrealistic  expectations  on  the 
system. 

A  first  consideration  in  setting  up  a  sediment  budget  is  to  define  the 
quantities  desired  (i.e.  unknowns)  and  the  quantities  that  are  assumed 
known  (i.e.  measured  or  calculated).  To  establish  the  procedure  fol¬ 
lowed  in  this  technical  note,  a  simplified  first  example  will  be  provided 
as  shown  in  Fig.  1  which  is  a  3  cell  sediment  budget  with  4  assumed 
knowns  (i.e.  the  volume  change  within  the  cell  and  the  net  sediment 
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transport  at  the  downdrift  boundary  of  the  downdrift  cell)  and  3  un¬ 
knowns  (i.e.  entering  and  exiting  quantities  at  the  intermediate  bound¬ 
aries  of  the  cells  along  with  the  net  sediment  transport  at  the  updrift 
boundary  of  the  updrift  cell).  In  Fig.  1,  the  updrift  side  of  the  inlet  is 
represented  by  the  cell  on  the  left  (i.e.  cell  1),  while  the  downdrift 
side  of  the  inlet  is  represented  by  the  cell  on  the  right  (i.e.  cell  3).  The 
middle  cell  (cell  2)  is  the  cell  with  the  inlet  in  it.  In  this  first  example, 
using  one  equation  for  each  cell,  continuity  balance  allows  for  the  3 


equations  as  follows: 

Qnu-Qsu  = ^ 

(la) 

Qsu-Qsd 

(lb) 

Q-sd  Q-nd 

(lc) 

where  Qy  represents  the  sediment  transport  in/out  of  the  cell  with 
the  subscripts  “n”  referring  to  net,  “u”  referring  to  updrift,  “s’  referring 
to  shoal/inlet,  “d”  referring  to  downdrift,  and  where  ^  represents  the 
volume  change  within  the  "i"  cell. 

As  might  be  expected,  since  we  have  3  equations  and  3  unknowns, 
this  system  can  be  solved  to  provide  an  “exact"  answer  (albeit  based 
on  the  assumptions  made).  In  the  present  note,  the  equations  are  set 
up  in  matrix  form,  the  reasons  for  which  will  become  clear  in  further 
“seemingly  simple"  examples.  Upon  establishing  the  matrix  form  with 
the  “unknowns’  on  the  left  side  of  the  equations  and  the  “knowns”  on 
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Updrift  Cell  #1  Inlet  Cell  #2  Downdrift  Cell  #3 

Fig.  1.  Three  cell  sediment  budget  without  bypassing. 


the  right  side  of  the  equations,  the  matrix  system  can  be  shown  to  be  of 
the  following  form: 


or  in  the  same  matrix  shorthand  form  as  before  (although  with  mod¬ 
ified  matrix  definitions): 


1-1  O' 

0  1  -1 

<Zu 

0  0  1 

.  Qw . 

dV, 

dt 

dV2 

dt 

d-h  +  Q 

dt  y"' 


(2) 


and  using  matrix  shorthand  notation  for  the  above  quantities  the 
equation  system  can  be  written  as: 


A  - Q  =  RHS 


In  this  situation  as  the  equation  system  is  overdetermined  (i.e. 
more  equations  than  unknowns)  and  correspondingly  the  A  matrix 
is  not  square,  so  the  “best”  solution  is  a  “least  squares”  solution  pro¬ 
vided  by: 

Q=  (at-a)_1-(at-KHS)  (8) 


A  ■ Q  =  RHS  (3) 

where  A  is  the  matrix  of  0’s,  l's,  and  —  l's,  Q  is  the  vector  of  “un¬ 
knowns”,  and  RHS  is  the  right  hand  side  matrix  of  “knowns”.  The 
solution  (for  the  “unknowns")  to  this  expression  can  be  simply  stated 
as: 

Q_=A~'-RHS  (4) 


where  AT  is  the  matrix  transpose  of  A  and  where  it  is  assumed  that  the 
inverse  of  (AT -A)  exists  (which  it  does  in  this  case  and  can  be  found 
as 


iV'-A')  1  = 


2  1 

3  3 
1  2 

3  3  J 


where  A“ 1  is  the  inverse  of  A.  As  long  as  A  has  an  inverse,  then  the 
equations  can  be  solved  exactly  and  in  this  case  the  inverse  of  A  can 
be  easily  found  as: 


A-1 


'i  i  r 
0  1  1 
o  o  1 


(5) 


At  this  point  it  is  tempting  to  imagine  that  if  we  could  add  an  ad¬ 
ditional  equation  using  the  overall  sediment  budget  (of  all  3  cells) 
then  we  would  have  4  equations  and  4  unknowns  and  thus  be  able 
to  solve  all  the  unknowns  in  our  problem.  The  overall  constraint 
equation  for  the  three  cells  combined  can  be  found  as: 


Qnu  —  Qnd 


dVj  dV2  dV3 
dt  +  dt  +  dt 


(9) 


Let's  assume  now  that  for  the  same  sediment  budget  system,  that 
both  the  net  updrift  sediment  transport  Q„w  and  the  net  downdrift  sed¬ 
iment  transport  Qna  are  known.  We  still  have  three  equations  for  the 
three  sediment  cells  but  we  now  only  have  2  unknowns.  This  presents 
us  with  a  possibility  of  dropping  one  of  the  equations  and  getting  an 
answer  with  only  2  equations  and  2  unknowns,  but  an  advantage  of 
keeping  all  three  equations  can  be  gained  via  using  a  least  squares 
solution,  therefore  allowing  some  possible  error  in  the  quantities  as¬ 
sumed,  but  yet  giving  us  a  “best”  answer  to  our  problem.  For  this  case 
the  following  equations  are  used: 


-^_n 

dt  Unu 

(6a) 

o  - dv 2 

Qsd  _  dt 

(6b) 

dV3  +  o 
w+Q-nd 

(6c) 

where  again,  the  “unknowns”  are  on  the  left  side  and  “knowns”  on 
the  right  side  of  the  equations.  This  set  of  equations  can  be  written 
in  matrix  form  as: 


-1  O' 

1  -1 

0  1 

.  Qw . 

dVj 


dt 


-Q-n, 


-  dt 


dV2 

dt 

-  +  Q-nd 


(7) 


which  leads  to  a  matrix  expression  to  be  solved  as: 


r  1  -1  0  0 

Q-UU 

0  1  -10 

Qsu 

0  0  1-1 

Q.sd 

10  0-1 

-  Q-nd  - 

dV, 

dt 

dV2 

dt 

dV3 

dt 

dVj  dV 2 

- -  4 - -  + 

dt  dt 


dV3 
dt  . 


(10) 


or  using  matrix  shorthand  notation  as  before  (although  with  modi¬ 
fied  matrix  definitions): 


A  ■  Q  =  RHS 


In  this  case  the  inverse  of  A  is  singular  though  as  might  be  expected 
since  the  equations  are  not  independent!  As  is  readily  seen,  the  4th 
equation  is  a  simple  combination  of  the  first  three  equations  (i.e.  add¬ 
ing  the  first  3  equations  leads  to  the  4th  equation  and  no  information 
is  gained).  The  point  to  be  made  is  that  in  sediment  budgets  of  3  cells 
as  shown  here,  one  cannot  expect  to  solve  for  all  unknowns  and  added 
information  must  be  provided  (perhaps  the  ratio  of  net  drifts  at  the 
left  and  right  ends  of  the  system). 

At  this  point  it  is  instructive  to  look  at  another  aspect  of  sediment 
budgets  that  is  important  to  the  subject  of  natural  sediment  bypassing 
on  the  outer  bar  of  an  inlet  and  it's  implications.  This  situation  is  pre¬ 
sented  in  Fig.  2  which  is  similar  to  Fig.  1  but  with  the  assumption  of  a 
natural  bypassing  route  for  sediment  on  the  outer  ebb  shoal  of  the 
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Updrift  Cell  #1  Inlet  Cell  #2  Downdrift  Cell  #3 

Fig.  2.  Three  cell  sediment  budget  with  bypassing. 


inlet.  In  this  case  the  sediment  rate  assumed  to  be  directly  bypassed  on 
the  outer  bar  (from  the  updrift  to  the  downdrift  beaches)  is  given  by 
Qbp  and  where  other  quantities  are  provided  as  before.  Both  updrift 
and  downdrift  net  sediment  transport  quantities  are  known  at  the 
ends  of  the  system  and  the  intent  is  to  find  the  3  unknowns  provided 
on  the  left  hand  side  of  the  matrix  equation  given  below: 


1  0  1 

Q™' 

0  1  1 

Q*i 

1 

O 

.  Qi>p  . 

Q  -d-Yl 

CL„„  dt 


dV3 

dt 


+  Q™ 


(ii) 


where  the  lines  of  the  respective  matrix  equation  represent  the  con¬ 
tinuity  balance  of  the  updrift,  downdrift,  and  inlet  cells  respectively. 
Again  the  system  can  be  expressed  in  shorthand  notation  as: 


A  ■  Q  =  RHS 


quantities  in  the  sediment  budget.  Additional  assessment  has  been 
made  showing  where  equations  can  be  solved  and  where  added 
assumptions  or  constraints  need  to  be  made  (or  clarified)  in  order  to 
solve  a  sediment  budget  at  an  inlet.  Specific  informational  constraints 
regarding  natural  bypassing  and  net  transports  at  boundaries  have 
been  addressed  and  some  limitations  on  what  information  can  be 
derived  from  a  sediment  budget  addressed.  It  is  hoped  that  the  value 
found  in  the  approach  utilized  herein  will  encourage  both  engineers 
and  geologists  to  clearly  state  sediment  budget  assumptions  and 
correctly  address  sediment  budgets  such  that  best  management  of 
sand  resources  will  be  practiced  at  inlets. 
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Attempting  to  solve  this  equation  again  leads  to  the  fact  that  the  A 
matrix  is  singular,  therefore  the  solution  cannot  be  found  without  fur¬ 
ther  constraints  such  as  either  Qsu  or  Qsd  are  zero.  The  point  to  be 
made  is  that  even  such  a  simple  system  as  provided  above  cannot  be 
solved  to  assess  the  inlet’s  natural  bypassing  without  careful  assess¬ 
ment  of  required  assumptions. 

2.  Summary 

Methods  have  been  presented  to  assess  the  possibility  of  solving 
sediment  budget  unknowns  via  sediment  budget  matrix  equations 
without  prior  knowledge  of  magnitudes  or  directions  of  sediment  rate 
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